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Abstract. We investigate whether the observed high number of carbon- and nitrogen-enhanced extremely metal- 
poor stars could be explained by peculiar evolutionary properties during the core He flash at the tip of the red 
giant branch. For this purpose we compute a series of detailed stellar models expanding upon our previous work; 
in particular, we investigate if during the major He flash the penetration of the helium convective zone into 
the overlying hydrogen-rich layers can produce carbon- and nitrogen-rich abundances in agreement with current 
spectroscopic observations. The dependence of this phenomenon on selected model input parameters, such as 
initial metallicity and treatment of convection is examined in detail. 
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1. Introduction 

Spectroscopic observations of ultra metal-poor (UMP) 
stars have disclosed that carbon and nitrogen are signifi- 
cantly overabundant with respect to Fe, compared to typ- 
ic al halo stars w hich show [C/Fe]R:![N/F e]~0.0. According 
to [Beers (1999"| ) and [Rossi et al. (1999[ ) about 10% of the 



stars with [Fe/H]<-2.5 display [C/Fe]>1.0; this fraction 
rises to about 25% when [Fe/H]<— 3.0. Also N appears to 
be overabundant by the same amount, while a-elements 
show a ratio with respect to Fe that is typical of halo 
stars, i.e. [a/Fe]~0.4. As discussed, e.g., by Norris ct al 



(1997a), it is difficult to explain these abundance ratios 



in terms of binary star evolution or ejecta from very mas 
sive Population III supernovae, which in particular pro 
duce few nitrogen but vast amounts of oxygen (Heger & 
Woosl ev. 2002 ). so that the origin of the surface abun- 



dances l of UMP stars is still shrouded in mystery. 

Concerning the nature of the UMPs, different ideas 
may be followed, the first one being that they are true 
Pop III stars — i.e. of initially zero metallicity — with the 
observed heavy elements (Fe, etc.) resulting from pollution 
of just the envelope by other stars such as supernovae. In 
two recent papers ( Weiss et al. 2000| , Paper I, and Schlattl 



et al. 2001, Paper II) we have thus discussed in detail the 



evolution of low-mass zero-metal stars, including the effect 



of atomic diffusion, mass loss and surface metal pollution. 
One important result was that the small entropy barrier 
between the H- and He-rich regions allows the He-flash 
driven convective zone to penetrate the overlying H-rich 
layers at the tip of the red giant branch. The consequent 
inward migration of protons into high-temperature regions 
leads to an H-shell flash, which further increases the ex- 
tension of the central convective region. At the late stages 
of this phase the convective envelope deepens and merges 
with the inner convective zone. The initially metal-free 
surface is therefore enriched by a large amount of matter 
processed in He- and H-burning reactions. In particular, 
very high C and N abundances follow. We denote this 
process of surface C and N enrichment as HElium Flash 
induced Mixing (HEFM). 



This scenario — already suggested by Fujimoto et al 
(2000D — could potentially explain the anomalous abun 
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dance pattern at the surface of UMP stars, since via this 
mechanism high abundances of C and N are produced 
in the post He-flash phases. In Paper II we have shown 
that an initially metal-free star of 0.82 Af©, polluted by 
0.0003 Mq of Z=0.02 material, undergoes HEFM and re- 
produces in its post He-flash phase approximately lumi- 
nosity and surface gravity of two of the best studied UMP 
stars, namely CS 22892-052 and CS 22957-027. The stel- 
lar mass was chosen such that the age on the RGB of 
13.7 Gyr is compatible with the current estimates of the 
age of the universe, and with the age of CS 22892-052 es- 
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Table 1. Data for the UMP stars discussed in the 
gravity, effective temperature, [Fe/H], [C/Fe], [N/Fe], 
of CS 22892-052 have been estimated by us. 



text. The columns show star name, logarithm of the surface 
and ^^C/^'^C ratio, respectively. The errors on T^e and log(g) 



star 



log(g) 



Fo/H] 



[C/Fc] 



[N/Fc] 



CS 22892-055 


a 


1.5±0.5 


4850±100 


-2.97±0.20 


1.10±0.23 


1.0±0.52 




CS 22957-027 




2.25±1.0 


4839±130 


-3.43±0.12 


2.20±0.30 


2.0±0.50 




CS 22948-027 


c 


1.0±0.3 


4600±100 


-2.57±0.23 


2.0±0.18 


1.8±0.24 


«10 


CS 22949-037 


d 


1.7±0.3 


4900±100 


-3.79±0.16 


1.05±0.20 


2.7±0.40 





" data from Norris ct al. (1997a ) 

' C & N ab undances from jSIorris et al. 
data from |Aoki et al. (2002|) 



(1997b| ), [Fe/H], g, and T^h from [Bonifacio et al. (l"99^ 



data from Nforris et al. (2001) 



timated by Cowan et al. (1997 ) from nuclear chronology 
(15.6±4.6 Gyr). The amoimt and composition of the pol- 
luting material was adjusted such that the observed sur- 
face [Fe/H] abundance ([Fe/H]w— 3) was obtained and, as 
the result of the HEFM, [C/Fe] and [N/Fe] ratios consid- 
erably larger than zero resulted. However, the values of 
[C/Fe] and [N/Fe] predicted by the models are about 2 
orders of magnitude higher than the observed ones; i.e., 
we predict [C/Fe] w [N/Fe] «4, while observations of these 
stars yield [C/Fe] w [N/Fe] « 1-2 (see Table |l|). We note, 
however, that higher N overabundances are observed, e.g., 
in CS 22949-037 ([Norris et al., 20011). 



An alternative explanation for the C and N enhance- 
ment in UMPs is provided by HEFM during the ther- 
mally pulsating asymptotic giant-branch phase (TP-AGB) 
which occurs in metal-free stars with M^l M© (Cassisi 



ion and the accretion of further material, e.g., a nearby 
supernova. 

In summary, HEFM during the major He flash is 
presently the most promising astrophysical explanation 
for the high [C,N/Fe] ratios observed in UMP stars de- 
spite their failure of reproducing the absolute value of the 
C and N enrichment. Therefore we are studying to what 
extend HEFM occurs also in low-mass stars with initial 
Z>0, and whether those objects are showing the observed 
[C,N/Fe] ratios of CS 22892-052 and CS 22957-027. This 
extension of our model assumption is warranted by the 
results about star formation in primordial environments, 
which agree on the fact that the formation of (very) mas- 
sive objects is strongly favoured in the absence of metals 



(Abel et al., 2002; Bromm et al., 2002). If this indeed is 



et al., L996; Fujimoto et al., 2000). However, age, surface 



gr avit} |, and offoctivo temperature of thooo AGS otaro arc 



the case, the UMPs would constitute the extreme end of 
Pop II, and had formed out of material alrea dy bearing 
the signature of ind ividual supernova events ( Bhigeyama 



more d ifficult to reconcile with thooo of the oboorvod ob 
jects, which are therefore unlikely to be TP-AGB stars 
themselves. The idea in this case is that the observed 
low-mass UMP objects form binary systems with more 
massive AGB stars, which have transferred par t of their 
C- and N-e nriched envelope to their companion ( Fujimotc 
et al., ^ 000 ). A potential problem of this hypothesis is that 
HEFM during the TP-AGB of intermediate-mass stars 
may enhance the surface oxygen conten t similar to carbon . 
This feature can be found in models of ^iess et al. (2002| ) , 
which however do not agree in this respect with computa- 
tions of Chieffi et al. (2001 ). An oxygen enhancement of 
[0/Fe]wl-2 might not necessarily be in contradiction to 
observations, as only very weak upper limits on the [ O /Fe] 
ratio exist (Sneden et al., 1996; Norris et al., 1997a). But 



Tsujimoto, 1998[ ) 



the main problem of this scenario is that not all UMPs 
are in binary systems, and the orbital periods of those who 
are, are inconsistent with an AGB mass transfer paradigm 
(Preston & Sneden, 2001). Moreover, the origin of Fe in 
these stars is as uncertain as in our favoured scenario. It is 
presently not clear whether intermediate-mass stars with 
[Fe/H]«-3 undergo HEFM during the TP-AGB or rather 
evolve like ordinary Pop II objects. In the latter case the 
surface C, N, O, and Fe abundances of UMP stars would 
then be the result of mass transfer from their compan- 



The final alternative for the origin of the carbon-rich 
UMPs is that they were already born with the abundances 
they now show. This scenario would receive support if 
single dwarfs of similar composition should definitely be 
identified. While one shifts the problem of identifying the 
source of the peculiar abundances to other stars, it is 
worthwhile, nevertheless, to investigate the evolution of 
such C-rich stars. For example, they are expected to reach 
standard RGB-tip luminosities and low gravities, which 
could be compared with those of observed objects. Stars 
of extremely low metallicity (Z = in particular), have a 
definitely shorter RGB, in contrast. 

In addition to the three plausible models for the ori- 
gin of the C-rich UMPs, we also investigate in this paper 
a physical problem associated with the HEFM. Due to 
the large uncertainties still affecting the treatment of con- 
vection in stellar models, we will assess the sensitivity of 
the HEFM process to different assumptions about the ex- 
tension and mixing efficiency of convective regions in the 
stellar interiors. 

The paper is structured as follows. Section 2 sum- 
marizes briefly the model computational techniques, in- 
put physics, and evolution of low-mass stars undergoing 
HEFM. Section 3 discusses the occurrence of HEFM and 
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Fig. 1. Evolution in the gravity- 
effective temperature diagram of mod- 
els with different initial metallicities: a 
M=0.82 M0-models with initial Z=0 
and surface metal pollution (dash- 
dot-dot-dotted line), a 0.78 Af0-model 
with Z=2.0 X 10"^ ('best model', 
solid line), and a O.82M0-model with 
Z=6.0 X 10"'*, where an initial C-, N- 
enhanced metal-mixture has been as- 
sumed (dashed line). In addition, a 
model with Z^2.0 x 10"^ and addi- 
tional overshooting {F = 0.03) un- 
dergoing a HEFM during the TP- 
AGB phase is shown (dash-dotted line) . 
Open symbols show the position of the 
two UMP stars CS 22892-052 and CS 
22957-027 (Table |^. The vertical error 
bars are slightly offset in order to not 
overlap. 



resulting surface abundances for different initial metallici- 
ties and various parameterizations of the stellar convective 



2 higher than the rates of Caughlan fc Fowler (198g ) and 



regions; Sect. 4 summarizes the results. 



about 30-50 % higher than the most recent value of 
et al. (2002| ). 



Kun2 



2. Modelling and evolution of stars undergoing 
HEFM during the He flash 

All calculations (as the ones in Paper II) have been per- 



formed using the Garching stellar evolution code (Weiss 
fc Schlattl 2000| ). The numerical details are presented in 



Paper II and will not be repeated here. The important fea- 
ture of our code is the ability to follow in detail the evolu- 
tion through the He flash considering simultaneously mix- 
ing and burning processes. Mixing is described by means 
of a time-dependent algorithm that treats convection as 
a fast diffusive process, where the diffusion constant is 
proportional to the convective velocity obtained from a 
convection theory. We have discussed in Paper II that the 
convective velocities predicted by both mixing-length and 



Canuto fc Mazzitelli's (1991) theory lead to very similar 



results. In general, arbitrary changes by up to 2 orders of 
magnitude in the mixing speed do not affect appreciably 
the HEFM process. 



Atomic diffusion is tre ated according to Thoul et al 
(1994), opacities are from [glcsias fc Rogers (199(]| ) and 



Alexander fc Ferguson (1994), neutrino energy losses from 



Munakata et al. (198E), and the equation of state (EOS) 



is a Saha EOS plus a simplified description of the degener- 



ate electron gas in the core regions following Kippenhahn 



fc Weigert (199C). Our reaction network follows the abun- 



dances ofH, ^He, ^He, ^^C, ^^C, ^^N, ^^N, i^o, i^q, with- 
out assuming a priori equilibrium compositions for any of 
these chemical elements. The reaction rates for the pp- 



and CNO-burning have been taken from Adelberger et al 
(1998), while the ^^C(q!, 7)^^0-rate agrees with Caughlan 



et al.'s (1985) value. The latter one is, in the relevant tem- 
perature range (0.1 < r/(10^ K) < 0.2), about a factor of 



Figure]^, as a reference, displays the evolutionary track 
of the model with 0.82 M©, Z=0, and surface chemical 
pollution mentioned in the introduction (similar to model 
Bl of Paper II), which experiences the HEFM. 

The 0.0003 Mq of Z=0.02 polluting material are in- 
stantaneously accreted before the star reaches the zero- 
age main sequence; this approach maximizes the effect of 
accretion and diffusion effects. After this episode the star 
evolves along the main sequence producing energy through 
the pp-chain. Diffusion is never able to bring the surface 
metals deep enough to contribute to the burning. 

At the sub-giant branch the core increases its temper- 
ature and density, and gradually 3-a reactions are setting 
in. In stars more massive than w 0.9 M0 the central H 
content is not fully exhausted when sufficient C is pro- 
duced (X(^^C) « 10~^°) to gain considerable nuclear en- 
ergy from the CNO-cycle. As a consequence a thermal 
runaway through the CNO-cycle ensues, causing a char- 
acteristic blue loop in the H-R diagram. The consequent 
expansion of the inner regions reduces the efficiency of 
the 3-a process and CNO-cycle, terminating the runaway 
which lasts only about 10^ years. In less massive models, 
as the ones considered in this work, the central H is al- 
ready exhausted before sufficient C can be produced, and 
thus no blue-loop occurs on the sub-giant branch. 

Regardless of the occurrence of this CNO flash, all low- 
mass models settle on the RGB and evolve towards higher 
luminosities. While on the sub-giant branch their main 
energy source remains the pp-chain, the contribution of 
the CNO-cycle is increasing gradually with luminosity as 
more carbon is created at the inner tail of the H-burning 
shell by 3-a reactions. In a 0.8 Mq star about 50% of the 
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H-burning energy star is produced by the CNO-cycle at 
the tip of the RGB. 

At this stage the mass difference between the location 
of the maximum energy generation in the core and the 
He-core boundary is Mq. Soon after the onset of 

the He flash a convective shell develops, which reaches the 
H-rich matter in the envelope about 1 month after the 
flash (left arrow in Fig. ||); as soon as H is carried into 
the interior hotter convective region, it starts to burn at a 
very high rate. Due to this extra energy input the upper 
boundary of the convective shell immediately moves closer 
to the surface (thus ingesting even more protons), while 
the He-burning rate is significantly reduced. The single 
inner convective region splits into two zones, one for the 
He- and one for the H-burning. The rapidly weakening He- 
burning shell can no longer support the underlying layers 
against contraction, and the released gravothermal energy 
forms a further convective region which disappears again 
after a few 1,000 year^ 

About 500 yr after the H ingestion the convective enve- 
lope deepens and merges with the convective region above 
the H-burning region (horizontal arrow in Fig. y). A huge 
amount of matter processed through H and He burning is 
brought to the surface enriched in He, C (produced during 
the He burning) and N (produced by the CNO cycle at the 
expenses of C). The resulting surface [C/Fe] and [N/Fe] 
ratios are, respectively, 4.1 and 4.3, while the surface oxy- 
gen abundance is practically unaltered with respect to the 
initial value (Table ^ . 

As a result of this dredge up of heavy elements (which 
happens on timescales of a few weeks) the envelope opacity 
is increased significantly; this produces an abrupt discon- 
tinuity in the effective temperature visible in Fig. ^ At 
this stage, with a practically ceased He burning, the star 
behaves as a newborn RGB star, climbing its own RGB. 
The H shell produces the energy needed to support the 
star. Another He flash ensues at the end of this second 
RGB phase, this time with a smaller and less degener- 
ate core. The flash is weaker and no HEFM happens. The 
subsequent evolution does not show any peculiar feature. 

As a general rule, all the changes in the physical in- 
puts and/or initial conditions that can contribute to an 
increase of the electron degeneracy and thus cause the lo- 
cation of the He-flash ignition to move closer to the border 
of the He core favour the occurrence of HEFM. Changes 
in the relative location as small as «0.02 Mq can make 
a difference between models undergoing HEFM and mod- 
els which do not experience it. In Paper II we have shown 
how an increase of the initial He abundance, or an increase 
of the initial mass, or the inclusion of heavy element pol- 
lution at the stellar surface are disfavouring the onset of 
HEFM. On the other hand, the inclusion of element diffu- 
sion favours this process, while the inclusion of mass loss 



^ This convective region appeared also in the calculations 
presented in Paper 11, but was not shown in the corresponding 
figure due do insufficient plot data. 
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Fig. 2. Evolution of the different convective zones — in- 
dicated by shaded areas — from the He-flash ignition on- 
wards in the Z=0 model of Fig. |l| with surface pollution. 
The dashed and solid lines represent the mass shells (in 
units of fraction of total mass) of maximum energy re- 
lease by H and He burning, respectively. The two arrows 
mark the penetration of the He-flash driven convective 
zone into H-rich layers and the merger of the convective 
envelope with the C- and N-rich inner convective region, 
respectively. The moment of maximum He-fusion energy 
release is indicated by the vertical dotted line. The zero 
point in time has been chosen to be 1 day before this 
moment. The short- living small radiative layers between 
the shells are only marginally stable against convection 
( "semi-convection" ) . 



from the stellar surface does not affect at all the onset of 
the HEFM. 

3. The influence of initial metallicity and 
convection 

In Paper II we have investigated the dependence of the 
HEFM process on several input-physics parameters and 
various assumptions adopted in computing stellar models; 
however, we have considered only stars with initial zero 
metallicity, plus some eventual metal pollution at the sur- 
face. It is worthwhile to investigate the dependence of the 
HEFM on the value of the initial stellar metallicity, given 
the unknown composition of these stars at the time of 
formation, as outlined in the introduction. 

Because of the uncertainty in the treatment of convec- 
tion in stellar structures, it is also important to verify how 
much the HEFM process is sensitive to different assump- 
tions about the extension and mixing efficiency of convec- 
tive regions in the stellar interior. In particular, we consid- 
ered various amounts of overshooting from the canonical 
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Fig. 3. Mass location of the He-core upper boundary and 
energy generation maximum as a function of the time be- 
fore the flash {to indicates the onset of the flash) , for the 
same models as in Fig. 0, indicated by the same line styles. 



formal convective boundaries flxed by the Schwarzschild 
criterion. 

In t hic eoction wo will diecuee thoeo offocte in thooroti — 



cal mo dole roproeontativo of CS 23802 053 and CS 22957 
027. 



3.1. Models with Z>0 

Models with initial y=0.23 and increasing Z have been 
computed in order to derive the value of the maximum 
initial metallicity which produces HEFM in a star with 
0.82 Mq. The heavy element mixture has been considered 
to be a-enhanced, with [Q;/Fe]=:0.4; carbon and nitrogen 
were scaled solar. We obtain HEFM only for Z < 10~^, 
corresponding to [Fe/H]<— 5.6. Since all observed [Fe/H]- 
abundances in UMPs so far are higher than —4, this 
already excludes the possibility that the C/N-anomalies 
were produced in stars with initial Fe-abundances as ob- 
served today. Due to the mixing of H-rich material into 
the H-depleted interior during the HEFM, the surface H 
abundance drops by about a factor of two, and the up- 
per limit for the surface Fe abundance after this phase is 
[Fe/H]<— 5.3, about 2 orders of magnitude smaller than 
the observed metallicities of CS 22892-052 and CS 22957- 
027. Whenever HEFM happens, the [C/Fe] and [N/Fe] 
ratios in the post He-flash phase are always of the order 
of 4, still about 2 orders of magnitude higher than the 
observed values. 

This upper value for the metallicity is slightly in- 
creased if we take into account the revised plasma- 



after the HEFM, i.e., 1 dex higher than in case of the 
old neutrino losses. The reason is that the degeneracy of 
the core increases and the He flash starts more off-centre. 
We also wish to notice here that, in spite of their differ- 
ent input physics and cruder treatment of the coupling 
between convection and burning, Fujimoto et al. (2000[ ) 
found a similar upper limit for the metallicity of a 0.80 Mq 
star undergoing HEFM, and similar values of the surface 
[C,N/Fe] ratios. 

We have also tes ted the effect of implementing a new 
EOSby Fwin (2000| ), which covers the entire stellar struc- 
ture in all relevant evol utionary phases, an d closely repro- 
duces the OPAL EOS ( [Rogers et al., 1996|) in the common 
validity range (see, e.g., ^chlattl, 2002| ). The inclusion of 
this new EOS in the models does not modify the value of 
the largest critical heavy elements abundance for which 
the HEFM occurs. However, with this equation of state 
the stellar ages are slightly reduced (by about 0.5 Gyr) 
with respect to our reference EOS, and therefore one has 
to choose a slightly lower mass to reproduce the observa- 
tions, which might favour the occurrence of HEFM. 

Therefore, we have computed a, what we call, 'best 
model' with an initial [Fe/H] equal to the observed one, 
and a mass which yields an age at the RGB tip well match- 
ing the age estimated for CS 22892-052 and CS 22957- 
027; we considered M=O.78M0, Z=2.0 x 10"^ (which 
corresponds, including the a enhancement by 0.4 dex, 
to [Fe / H] =— 3 .3) , neu trino emission following Haft et al 
(1994|) , llrwinl 's (|2000| ) EOS, plus atomic diffusion of H 



He, and metals (which further reduces the ages by about 
0.5 Gyr). The corresponding evolutionary track is shown 
in Fig. 1^, too (solid line). In spite of the reduced stellar 
mass the effect of the initial metallicity dominates, and no 
HEFM occurs. The location of the maximum energy gen- 
eration at the flash is 0.30 away from the boundary of 
the He core (solid line in Fig. y), a distance about 0.01 Mq 
larger than in the cases when HEFM occurs. We computed 
also the evolution of stars of 0.75 Mq and 0.82 Mq with 
the same initial composition as the 0.78 Mq star to span 
the age range determined for CS 22892-052 from nuclear 
chronology, not obtaining HEFM in any case. 



3.2. Models including overshooting 

Models including overshooting have been computed in or- 
der to check its influence on both the onset of HEFM 
and the amount of C and N dredged up to the surface. 
Since overshooting increases the extension of convective 
regions, the occurrence of HEFM is in principle favoured. 



neutrino rates by Haft et al. (1994). In this case the limit- 
ing metallicity is Z — lO^'', corresponding to [Fe/H]=— 4.3 



Our overshooting description follows Blocker et al. (1998) 
and it is modelled as an exponential diffusive process. The 
diffusion constant of the overshoot region decays expo- 
nentially outside the Schwarzschild convective boundary, 
starting from the value assumed by the convective dif- 
fusion coefficient at the convective boundary; the decay 
length is F times the pressure scale height at the con- 
vective boundary. Hydrodynamical simulations of shallow 
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convective envelopes in A stars by Freytag et al. (1996) 
predict f =0.25±0.05, while for the overshooting from the 
convective envelopes of AGB stars Herwig et al. (1997) 
took F = 0.02. 

We have computed a series of models (adopting the 
same input physics as in Paper II) with different values for 
F, considering M=0.82Mo, y=0.23 and Z=2.0 x 10"^ 
which corresponds to [Fe/H]=— 3.3, approximately the 
[Fe/H] value determined for the two UMP stars under 
scrutiny. In these models the simpler EOS (Saha-EOS 
completed with degenerate electron gas in the deep in- 
terior) and the neutrino losses of Munakata et al. (1985) 
have been used. 

We have first considered overshooting from all the 
convective boundaries, starting from the main-sequence 
phase. As a general result, we found that the inclusion 
of overshooting during the main sequence increases the 
lifetime of the small convective core which develops at 
the beginning of the main sequence, when the ^He abun- 
dance gradually reaches nuclear equilibrium in the centre; 
this, in turn, causes an increase of the main-sequence life- 
time and a decreased electron degeneracy in the He core 
along the RGB, which disfavours the occurrence of HEFM. 
Moreover, the increased lifetime would demand a larger 
stellar mass in order to obtain not too high stellar ages, 
and that would further reduce the probability for HEFM. 

Models with A/=0.82Mo and F=0.03 do not show 
HEFM, and are continuing their evolution through the 
horizontal branch to the TP-AGB (dash-dotted line in 
Fig. H). During this phase the inter-shell convective zone 
is able to penetrate the H-rich envelope, similar to the 
HEFM process ( [Cassisi et al., 1996| ; [Fujimoto et al., 2000| ). 



The surface C and N abundances rise again by about 4 dex, 
while O is increased by about 3 dex. Hence, an O enrich- 
ment about one order of magnitude smaller than for C is 
expected in these stars, which implies that [0/Fe] should 
be at most 1 for the stars in Table |l|. Since no oxygen 
abundance could be determined in these stars, but only 
an upper limit of [O/Fe]^0.6 (Sneden et al., 1996, for 



CS 22982-052), this scenario could provide a possible ex- 
planation for C-rich UMPs, too. However, surface gravity 
and temperature of these AGB stars are too low to be in 
agreement with the stars under scrutiny (Fig. 

Higher values of F do not lead to HEFM. In fact, for 
= 0.1 there is no He flash, but quiescent He-burning 
ignition and of course no HEFM. 

In case overshooting is included only from the RGB 
phase onwards, the stellar lifetime is not affected, and one 
needs at least F = 0.17 to get HEFM. This value is about 



30% smaller as the one obtained by Freytag et al. (1996) 
for A-star envelopes, but appears to be quite high com- 
pared to that needed in AGB stars (Herwig et al., 1997); 
moreover, switching on overshooting only during the RGB 
phase appears to be entirely ad hoc; in addition, the re- 
sulting [C/Fc] and [N/Fe] ratios after the HEFM are again 
of the order of 4, i.e., too high with respect to the obser- 
vations. 
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Fig. 4. Same model as in Fig. ^ but reducing the convec- 
tive mixing efhciency by a factor 2 x 10'*. 



3.3. Models with reduced mixing efficiency 

From the computations performed until now, it appears 
that one of the main problem to the HEFM scenario is the 
too high C and N enrichment of the surface. One solution 
to obtain lower surface [C,N/Fe] ratios would be to have 
a smaller overlap between the convective region develop- 
ing at the flash, and the overlying H-rich region. Smaller 
overlap implies less protons and therefore less processed C 
and N. This smaller overlap could be achieved by strongly 
reducing the diffusion coefficient associated to the convec- 
tive mixing. Using standard mixing-length velocities, the 
crossing time of the convective shell is of the order of a few 
hours. A reduction of the "convective" diffusion coefhcient 
is performed by simply reducing the multiplicative factor 
that links the coefficient itself to the convective velocity. 

The upper boundary of the He-flash driven convective 
shell is at first order dictated by the energy of the flash 
and the location of the ignition point. However, after some 
protons are engulfed the increased energy production (due 
to the additional H burning) pushes the upper boundary 
further up, thus increasing the number of proton ingested. 
If the mixing efficiency is reduced, protons are mixed less 
deep and less energy is produced, with a smaller increase of 
the convective shell extension. This would cause a smaller 
final number of protons ingested and a lower N production. 

In order to explore this scenario in more detail, we 
recomputed the polluted Z=0 model of Sect. 2, which 
experiences HEFM, decreasing the "convective" diffusion 
coefficient by a factor of 2 x 10"*. The extension of the con- 
vective regions are displayed in Fig. ^. The upper bound- 
ary of the H-flash driven convective shell extends less into 
H-rich layers, and the location of the H-burning region 
is less deep; the thickness of the H convective shell be- 
fore the dredge up is also thinner, implying a reduction 
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Table 2. Surface abundances of initial metal- free 
0.82 Mq stars polluted by 3 x 10"** Mq of alpha-enhanced 
([a/Fe]=0.4]) material with Z = 0.02 at the ZAMS. AU 



models contain updated neutrino- losses (Haft et al., 1994) 
and Irwin's ( |200C ) EOS. The evolution of the convective 
regions after the flash are shown in the figures denoted in 
the first column. Jd is the factor by which the standard 
diffusion efficiency has been multiplied in each model. 
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our tests in Sect. 2 there appear to be two main problems 
with the HEFM scenario: 

Firstly, too much of C and N is transported to the 
surface after the HEFM, and secondly, no HEFM appears 
to occur for an initial metallicity that can match the ob- 
served [Fe/H] of the stars under scrutiny, unless surface 
heavy element pollution has been efficient, as discussed in 
Paper II. 

Once again, we want to point out that, in spite of 
their different input physics and cruder treatment of 
the coupling between convection and burning, Fujimotc 



5x10" 



2.3 3.2 



0.4 



-3.4 5.6 0.07 



of the amount of C dredged up with respect to the case 
with standard convective mixing. The final surface abun- 
dance ratios [C/Fe] and [N/Fe] are in this case 2.3 and 3.2 
(Table H), respectively, closer to the observed values. 

The surface *^C/^'^C ratio is also affected, being now 
increased to 5.6 compared to 4.5 which we obtain by us- 
ing the standard mixing efficiency. This increased value is 
in slightly better agreement with the observed values of 
about 10. In spite of the much reduced mixing efficiency 
the dredge up still happens fast, on timescales of the order 
of lO^^ years. 

4. Discussion 

The results shown in the previous section highlight the 
fact that it is very difficult to reproduce the surface abun- 
dances of CS 22892-052 and CS 22957-027 from single-star 
evolution, at least in the case the observed surface abun- 
dances pattern is not primordial. To investigate this alter- 
native, we have computed additionally the evolution of an 
0.82 Mq star, with the same physics of our 'best model' 
and r=0.23, but with Z=& x 10"^ and [C/Fe] = [N/Fe] =2 
(implying [Fe/H]=— 3.3), in agreement with the observed 
abundances. The evolutionary track is shown in Fig. 0, 
while the He-core boundary before and at He ignition is 
plotted in Fig. |[ As expected, this model does not expe- 
rience HEFM (the He flash starts very deep in the core), 
so that the surface [C/Fe] and [N/Fe] ratios are not mod- 
ified with respect to the initial values. Note that the ef- 
fective temperature along the RGB is higher than for the 
case with HEFM, because of the much lower C-abundance. 
The location of the track fits reasonably well the observed 
gravities and effective temperatures, and thus, this model 
would be able to reproduce the observed properties of 
CS 22957-027. However, as discussed in Sect. 1, if this 
were a realistic model for the observed object, the source 
for the high carbon and nitrogen abundances still remains 
unknown, as it is not clear how to pollute the interstellar 
medium with matter showing these anomalous abundance 
ratios. 

Thus, assuming that the observed surface composition 
is not primordial, the HEFM is presently the only way 
to get a surface C and N enhancement without invoking 
some kind of mass transfer from binary companions. From 



et al. (2000) found similar upper limits for the metallic- 
ity of stars undergoing HEFM, and similar values of the 
[C,N/Fe] ratio for stellar mass and metallicity values sim- 
ilar to our 'best model'. 

Concerning the appearance of HEFM in stars with 
[Fe/H] in agreement with observation s, on e possible so- 
lution would be, as discussed in Sect. 3.3, a completely 
ad hoc choice of overshooting efficiency. Another possi- 
bility is related to the H depletion during the dredge up 
of the HEFM products. During this phase the convective 
envelope reaches deeper regions involved in the HEFM, 
which are devoid of H. According to our models the sur- 
face [Fe/H] increases by Ri0.3dex because of the dilution 
of H. If the drop in hydrogen abundance is about five times 
higher, a surface [Fe/H] of about —3.3 could be obtained 
with an initial [Fe/H]=— 4.3 corresponding to the upper 
limit of getting HEFM. 

In order to achieve a larger H depletion, a smaller enve- 
lope thickness (in mass) is needed, i.e., a smaller H reser- 
voir. We estimate that right before the flash ignition an 
envelope of only about 0.03 Mq is needed in order to re- 
produce the observed post-HEFM [Fe/H] surface values 
with our 'best model'; with this envelope mass the star 
would still be able to ignite He close to its RGB location 
(see, e.g., Castellani fc Castellani 1993| ), i.e., the effective 



temperature and gravity after the HEFM would still be in 
accordance with CS 22892-052 and CS 22957-027. Taking 
into account the age of the star of about 14Gyr, which 
fixes its initial mass to be about 0.8 Mq, it would imply 
that about 0.27 Mq are lost during the evolution, a value 
which is not much higher than the average mass loss ex- 
perienced by globular cluster stars (about 0.20 Mq). The 
lack of metals in the outer layers of UMPs, however, should 
cause a strong decrease of mass-loss efficiency in compar- 
ison with more metal-rich stars, if radiative driven winds 
are the main source of mass loss. Nevertheless, due to our 
poor knowledge about mass-loss processes in RGB stars, 
this possibility cannot be completely ruled out. 
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